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Abstract 



We present a systematic study of 400 combinations of the charged lepton and 
neutrino mass matrices with six vanishing entries or texture zeros. Only 24 
of them, which can be classified into a few distinct categories, are found to be 
compatible with current neutrino oscillation data at the 3a level. A peculiar 
feature of the lepton mass matrices in each category is that they have the 
same phenomenological consequences. Taking account of a simple seesaw 
scenario for six parallel patterns of the charged lepton and Dirac neutrino 
mass matrices with six zeros, we show that it is possible to fit the experimental 
data at or below the 2a level. In particular, the maximal atmospheric neutrino 
mixing can be reconciled with a strong neutrino mass hierarchy in the seesaw 
case. Numerical predictions are also obtained for the neutrino mass spectrum, 
flavor mixing angles, CP-violating phases and effective masses of the tritium 
beta decay and the neutrinoless double beta decay. 

PACS number(s): 12.15.Ff, 12.10.Kt 
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I. INTRODUCTION 



Recent solar [1], atmospheric [2], reactor (KamLAND [3] and CHOOZ [4]) and accelerator 
(K2K [5]) neutrino oscillation experiments have provided us with very convincing evidence 
that neutrinos are massive and lepton flavors are mixed. In the framework of three lepton 
families, a full description of the lepton mass spectra and flavor mixing at low energies needs 
twelve physical parameters: 

• three charged lepton masses m e , and m r , which have precisely been measured [6]; 

• three neutrino masses m 1? m 2 and m 3 , whose relative sizes (i.e., two independent 
mass-squared differences Am^ = m\ — m\ and Am| x = m\ — m 2 ) have roughly been 
known from solar (Am^ ~ 1CT 5 eV 2 ) and atmospheric ((Am^l ~ 1CT 3 eV 2 ) neutrino 
oscillations; 

• three flavor mixing angles #12, #23 and #13, whose values have been determined or 
constrained to an acceptable degree of accuracy from solar (# 12 ~ 33°), atmospheric 
(#23 ~ 45°) and reactor (9 13 < 13°) neutrino oscillations; 

• three CP-violating phases 5, p and a, which are completely unrestricted by current 
neutrino data. 

The future neutrino oscillation experiments are expected to fix the sign of Am^, to pin down 
the magnitude of #13 and to probe the "Dirac-type" CP-violating phase 5. The proposed 
precision experiments for the tritium beta decay and the neutrinoless double beta decay 
will help determine or constrain the absolute scale of three neutrino masses. Some informa- 
tion about the "Majorana-type" CP-violating phases p and a may also be achieved from a 
delicate measurement of the neutrinoless double beta decay. However, it seems hopeless to 
separately determine p and a from any conceivable sets of feasible neutrino experiments in 
the foreseeable future. 

The phenomenology of lepton masses and flavor mixing at low energies can be formulated 
in terms of the charged lepton mass matrix M; and the (effective) neutrino mass matrix M v . 
While the former is in general arbitrary, the latter must be a symmetric matrix required by 
the Majorana nature of three neutrino fields. Hence we diagonalize Mi by using two unitary 
matrices and M v by means of a single unitary matrix: 
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The lepton flavor mixing matrix V is defined as V = UjU u , which describes the mismatch 
between the diagonalizations of Mi and M v . In the flavor basis where Mi is diagonal and 
positive, V directly links the neutrino mass eigenstates (^1,^2,^3) to the neutrino flavor 
eigenstates (v e , u^, v T ): 
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A convenient parametrization of V is 
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sin 6^- (for = 12, 23, 13). We have known that #23 > 9 12 > 9 13 
holds, but how small 9 13 is remains an open question. A global analysis of current neutrino 
oscillation data shows that #13 is most likely to lie in the range 4° < 9± 3 < 6° [7] . In this case, 
we are left with a bi-large mixing pattern of V, which is quite different from the tri-small 
mixing pattern of the quark flavor mixing matrix. 

To interpret the observed hierarchy of Am^ and Ara^ as well as the bi-large lepton flavor 
mixing pattern, many phenomenological ansatze of lepton mass matrices have been proposed 
in the literature [8] . A very interesting category of the ansatze focus on the vanishing entries 
or texture zeros of Mi and M v in a specific flavor basis, from which some nontrivial and 
testable relations can be established between the flavor mixing parameters and the lepton 
mass ratios. We argue that texture zeros of lepton mass matrices might result from a kind 
of new (Abelian or non-Abelian) flavor symmetry beyond the standard electroweak model. 
Such zeros may dynamically mean that the corresponding matrix elements are sufficiently 
suppressed in comparison with their neighboring counterparts. From a phenomenological 
point of view, the study of possible texture zeros of M; and M v at low energies do make 
sense, because it ought to help reveal the underlying structures of leptonic Yukawa couplings 
at a superhigh energy scale. 

The main purpose of this article is to analyze the six- zero textures of Mi and M v in a 
systematic way. To be specific, we take Mi to be symmetric, just as M v is. This point is true 
in a number of SO (10) grand unification models, in which the group symmetry itself may 
dictate all fermion mass matrices to be symmetric [9]. Then a pair of off-diagonal texture 
zeros in Mi or M u can be counted as one zero. We further require that each mass matrix 
contain three texture zeros, such that the moduli of its three non-vanishing elements can 
fully be determined in terms of its three mass eigenvalues 1 . Because there exist 20 different 



One may certainly consider the possibility that one mass matrix contains two zeros and the other 
consists of four zeros. In this case, the former loses the calculability - namely, its four independent 
moduli cannot completely be calculated in terms of its three mass eigenvalues; and the latter causes 
the correlation between one of its three mass eigenvalues with the other two - this kind of mass 
correlation is in general incompatible with the relevant experimental data. One must reject the 
possibility that one mass matrix consists of one zero and the other contains five zeros, because the 
latter only has a single non-vanishing mass eigenvalue and is in strong conflict with our current 
knowledge about the charged lepton or neutrino masses. Therefore, we restrict ourselves to the 
most interesting and feasible case: six texture zeros are equally shared between M/ and M v . 
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patterns of Mi or M v with three texture zeros, we totally obtain 20 x 20 = 400 combinations 
of Mi and M v with six texture zeros. A careful analysis shows that only 24 of them, which 
can be classified into a few distinct categories, are consistent with current neutrino oscillation 
data at the 3a level. We find that the lepton mass matrices in each category have a peculiar 
feature: they do not have the same structures, but their phenomenological consequences 
are exactly the same. This isomeric character makes the six-zero textures of lepton mass 
matrices especially interesting for model building. It is noticed that those 24 patterns of 
Mi and M v are difficult to agree with today's experimental data at the 2a level, mainly due 
to a potential tension between the smallness of Ami^/IAmU and the largeness of sin 2 #23. 
Taking account of a very simple seesaw scenario for six parallel patterns of the charged lepton 
and Dirac neutrino mass matrices with six zeros, we demonstrate that it is possible to fit 
the present neutrino data at or below the 2a level. In particular, the maximal atmospheric 
neutrino mixing (i.e., sin 2 2# 23 ~ 1) can be reconciled with a strong neutrino mass hierarchy 
in the seesaw case. Specific numerical predictions are also obtained for the neutrino mass 
spectrum, flavor mixing angles, CP-violating phases and effective masses of the tritium beta 
decay and the neutrinoless double beta decay. 

The remaining part of this article is organized as follows. A classification of the six-zero 
textures of lepton mass matrices is presented in section II, where a few criteria to select 
the phenomenologically favorable patterns of Mi and M u are also outlined. Section III is 
devoted to the analytical and numerical calculations of 24 patterns of lepton mass matrices 
with or without the structural parallelism between Mi and M v . A simple application of the 
seesaw mechanism to the charged lepton and Dirac neutrino mass matrices with six texture 
zeros is illustrated in section IV. Finally, we summarize our main results in section V. 



II. A CLASSIFICATION OF THE SIX-ZERO TEXTURES 



A symmetric lepton mass matrix M (i.e., Mi or M u ) has six independent entries. If three 
of them are taken to be vanishing, we totally arrive at 
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patterns, which are structurally different from one another. These twenty patterns of M 
can be classified into four categories: 



1. Three diagonal matrix elements of M are all vanishing (type 0): 
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where those non- vanishing entries are simply symbolized by x's. 
2. Two diagonal matrix elements of M are vanishing (type I): 
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which are of rank three; and 
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which are of rank two. 
3. One diagonal matrix element of M is vanishing (type II): 
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which are of rank three; and 
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which are of rank two. 
4. Three diagonal matrix elements of M are all non- vanishing (type III): 



Mm = 



fx 0\ 

0x0 
v x j 



(9) 



(10) 



We see that Mo and Mi x are the well-known Zee [10] and Fritzsch [11] patterns of fermion 
mass matrices, respectively. Both of them are disfavored in the quark sector [12]. While the 
original Zee ansatz is also problematic in describing lepton masses and flavor mixing [13], 
the Fritzsch ansatz is found to be essentially compatible with current neutrino oscillation 
data [14]. 

Allowing the charged lepton or neutrino mass matrix to take one of the above three- zero 
textures, we totally have 20 x 20 = 400 combinations of Mi and M v . We find that 141 of 
them can easily be ruled out. First, the pattern in Eq. (5) is not suitable for M h because 
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three charged leptons have a strong mass hierarchy and the sum of their masses (i.e., the 
trace of M u ) cannot be zero. Second, the rank-two patterns in Eqs. (7) and (9) are not 
suitable for Mj, because the former must have one vanishing mass eigenvalue. Third, Mj 
and M y cannot simultaneously take the pattern in Eq. (10), otherwise there would be no 
lepton flavor mixing. We are therefore left with (20 — 7) x 20 — 1 = 259 combinations of Mj 
and M v . 

To pick out the phenomenologically favorable six- zero patterns of lepton mass matrices 
from 259 combinations of M; and M u , one has to confront their concrete predictions for the 
lepton mass spectra and flavor mixing angles with current neutrino oscillation data. The 
strategies to do so are outlined below: 

1. For each combination of Mi and M u , we do the diagonalization like Eq. (1). Because 
Mi has been specified to be symmetric, Ui = U* must hold. The matrix elements 
of Ui can be given in terms of two mass ratios (x l = m e /m^ m 0.00484 and y { = 
Trifi/rrir ~ 0.0594 [6]) and two irremovable phase parameters 2 . A similar treatment is 
applicable for the neutrino sector. The ratio of two independent neutrino mass-squared 
differences reads 



R v = 



Am 2 



21 



Am 2 



31 



I 1 • L vi)v\ 



where x v = mi/m 2 and y v = m 2 /m 3 . Note that x v < 1 (i.e., m\ < m 2 ) must hold, but 
it remains unclear whether y v < 1 (normal mass hierarchy) or y u > 1 (inverted mass 
hierarchy). The numerical results for Am^ and jAm^l, which are obtained from a 
global analysis of current neutrino oscillation data [15], have been listed in Table I. 
We are therefore able to figure out the allowed range of R v . 

2. The lepton flavor mixing matrix V = U\\J V can then be obtained. Its nine elements 
depend on four mass ratios {x x ,y L ,x u and y v ) and two irremovable phase combinations, 
which will subsequently be denoted as a and j3. In the standard parametrization of 
V, as shown in Eq. (2), one has 
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The experimental results for sin 2 12 , sin 2 9 2 3 and sin 2 13 are also listed in Table I. 
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Without loss of generality, we can always arrange one of the three non-vanishing entries of M; 
(or M v ) to be positive. We are then left with two free phase parameters in Mi (or M„). 
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3. With the help of current experimental data, we make use of Eqs. (11) and (12) to 
look for the parameter space of each pattern of lepton mass matrices. The relevant 
free parameters include two neutrino mass ratios [x v and y u ) and two CP-violating 
phases (a and (3). The latter may in general vary between and 2n. In our numerical 
analysis the points of x u , y u , a and (3 will be generated by scanning their possible 
ranges according to a flat random number distribution. Thus the density of output 
points in the (x u ,y u ) and (a,/3) plots will be a clear reflection of strong constraints, 
imposed by the neutrino oscillation data and the model (or ansatz) itself, on these 
parameters. A combination of M t and M v will be rejected, if its parameter space is 
found to be empty. 

Of course, whether the parameter space of a specific pattern of lepton mass matrices is 
empty or not depends on the confidence levels of relevant experimental data. We shall focus 
on the 2cr and 3<x intervals of Am^, Am^, sin 2 #i 2 , sin 2 6> 2 3 and sin 2 13 given in Ref. [15]. 
It is worth mentioning that a plain scan of the unknown parameters (x u ,y u ) and (a, (3) is 
empirically simple and conservative, provided the reasonable ranges of Am^ etc have been 
fixed. In this approximation the error bars of those observables need not be statistically 
treated. 

Examining all 259 combinations of the charged lepton and neutrino mass matrices is a 
lengthy but straightforward work. We find that only 24 of them, whose Mi and M v both 
belong to type I given in Eqs. (6) and (7), are compatible with current neutrino oscillation 
data at the 3<r level. The detailed analytical and numerical calculations of those 24 patterns 
will be presented in sections III and IV. 

Once the parameter space of a given pattern of lepton mass matrices is fixed, one may 
obtain some predictions for the neutrino mass spectrum and leptonic CP violation. For 
example, the absolute values of three neutrino masses can be determined as follows: 



Three CP-violating phases in the standard parametrization of V are also calculable. As for 
CP violation in neutrino-neutrino or antineutrino-antineutrino oscillations, its strength is 
measured by the Jarlskog invariant J [16]. The definition of J reads 



where the subscripts (a,b,c) and (i,j,k) run respectively over (e,/i,r) and (1,2,3). The 
magnitude of J depends on both (x„,y„) and (a, (3). If \J\ ~ 1% is achievable, then 
leptonic CP- and T- violating effects could be measured in a variety of long-baseline neutrino 
oscillation experiments [17] in the future. 
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In addition, interesting predictions can be achieved for the effective mass of the tritium 
beta decay (m) e and that of the neutrinoless double beta decay (m) ee : 



(m) 2 e = Y,( m *\ V ^\ 2 ) = ™l(xly 2 u \V el \ 2 + yl\V e2 \ 2 + \V e3 \ 2 ) , 



(m) ee = 



E (^) 
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x,yuV e \ + y u V e l + V e 2 3 



(15) 



The present experimental upper bound on (m) e is (m) e < 2.2 eV [6], while the sensitivity of 
the proposed KATRIN experiment is expected to reach (m) e ~ 0.3 eV [18]. In comparison, 
the upper limit {m) ee < 0.35 eV has been set by the Heidelberg-Moscow Collaboration [19] 
at the 90% confidence level 3 . The sensitivity of the next-generation experiments for the 
neutrinoless double beta decay is possible to reach (m) ee ~ 10 meV to 50 meV [21]. 



III. FAVORED PATTERNS OF LEPTON MASS MATRICES 

The 24 patterns of lepton mass matrices, which are found to be compatible with current 
neutrino oscillation data at the 3<r level, all belong to the type-I textures listed in Eqs. (6) 
and (7). To make our subsequent discussions more convenient and concrete, we rewrite those 
type-I textures of Mi or M v and list them in Table II. Two comments are in order. 

• Each type-I texture of M (i.e., Mi or M v ) can be decomposed into M = PMP T , 
where P denotes a diagonal phase matrix and M is a real mass matrix with three 
positive non-vanishing elements. The diagonalization of M requires an orthogonal 
transformation: 



0*MO* = A 2 , (16) 




where Aj (for i — 1, 2, 3) stand for the physical masses of charged leptons (i.e., Ai^ = 
m e,/i,r) or neutrinos (i.e., Aj = rrii). Then the unitary matrix U (i.e., Ui or U u ) used to 
diagonalize M takes the form U = PO. 

Note that the matrix elements of M and O can be determined in terms of Aj. This 
calculability allows us to express the rank-3 (or rank-2) patterns of M in a universal 
way, as shown in Table II. It turns out that the relation 

M In = E a M ll E^, (n = l,...,6) (17) 
holds for those rank-3 textures, where 



3 If the reported evidence for the existence of the neutrinoless double beta decay [20] is taken into 
account, one has 0.05 eV < (m) ee < 0.84 eV at the 95% confidence level. 
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As for three rank-2 textures, we have 

M h = /-.-..U, /•:.' , Mi 8 = E 4 M h Ej , 



M h = E 5 M lr E% . 



It is easy to check that E n is a real orthogonal matrix; i.e., E n E^ 
In addition, E A = E 2 E 3 = E 3 E 6 = E 6 E 2 and E 5 = Ej hold. 



(19) 
E 1 holds. 



Eqs. (17) and (19) will be useful to demonstrate the isomeric features of a few categories of 
lepton mass matrices with six texture zeros, as one can see later on. 



A. Six Parallel Patterns (rank-3) 

We have six parallel patterns of Mi and M u , 
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which are compatible with current neutrino oscillation data at the 3a level. Given M{" 
being diagonalized by the unitary matrix Ui tl/ , M^ (for n > 1) can then be diagonalized by 
E n Ui tU as a result of Eq. (17). The lepton flavor mixing matrix derived from M/ n and 
is found to be identical to V = UjU v , which is derived from M[ and Mf x : 

K = (E n Ui)H E M = Ul(E^E n )U u = V . (21) 

This simple relation implies that six parallel patterns of Mi and M v are isomeric - namely, 
they are structurally different from one another, but their predictions for lepton masses and 
flavor mixing are exactly the same [22]. It is therefore enough for us to consider only one 
of the six patterns in the subsequent discussions. With the help of Eq. (16), the moduli of 
three non-vanishing elements of M/ or M v are given by 

A = A 3 (l -y + xy) , 

\y(l-x)(l-y)(l+xy)] 1/2 

B = A3 1 , 

1-y + xy 

C = \J , Y\ (22) 

\l-y + xy) 

where the subscript "/" or "z/" has been omitted for simplicity. Furthermore, we obtain the 
matrix elements of O in terms of the mass ratios x and y (see Table II for the definition of 
\ and cj: 
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(23) 



Note that a 2 , b 2 and c 2 are imaginary, and their nontrivial phases arise from a minus sign of 
the determinant of M(i.e., Det(M) = — AC 2 e 2llp ). Because of < x v < 1 extracted from the 
solar neutrino oscillation data [1], we can obtain < y v < 1 from Eq. (22) as required by 
the positiveness of A v , B v and C u 4 . Hence the six isomeric patterns of lepton mass matrices 
under discussion guarantee a normal neutrino mass spectrum. 

Nine elements of the lepton flavor mixing matrix V 
be written as 



UIU V = Oj(PlP u )O u can explicitly 

V pq = (a l p )*a^ a + (b l p )*b^ + ( C J,)*< , (24) 

where the subscripts p and q run respectively over (e,/i, r) and (1,2,3), and the phase 
parameters a and (5 are defined by a = — <f t ) — /3 and /3 = (0^, — 4>i)- Note that F 
consists of four free parameters x u , y u , a and (3. The latter can be constrained, with the 
help of Eqs. (11) and (12), by using the experimental data listed in Table I (Am^ > as 
a consequences of < y u < 1). Once the parameter space of (x u ,y u ) and (a, (3) is fixed, 
one may quantitatively determine the Jarlskog invariant J and three CP-violating phases 
(5, p, a) . It is also possible to determine the neutrino mass spectrum and two effective 



4 Although y u > 1 is in principle allowed by rephasing the non- vanishing elements of M u , our 
numerical analysis indicates that this possibility is actually incompatible with current experimental 
data. 
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masses (m) e and (m) ee denned in Eq. (15). The results of our numerical calculations are 
summarized in Figs. 1-3. Some discussions are in order. 

1. We have noticed that the parameter space of (x u ,y u ) or (a, (5) will be empty, if the 
best-fit values or the 2a intervals of Am^, Am| 1; sin 2 #12, sin 2 623 and sin 2 #13 are 
taken into account. This situation is due to a potential conflict between the largeness 
of sin 2 #23 and the smallness of R u , which cannot simultaneously be fulfilled for six 
parallel patterns of M ; and M v at or below the 2a level. 

2. If the 3a intervals of Am^, Am^, sin 2 # 12 , sin 2 #23 an d sin 2 # 13 are used, however, the 
consequences of M t and M u on two neutrino mass-squared differences and three flavor 
mixing angles can be compatible with current experimental data. Fig. 1 shows the 
allowed parameter space of (x u ,y u ) and (a, (3) at the 3a level. We see that (3 ~ ir 
holds. This result is certainly consistent with the previous observation [14]. Because 
of y v ~ 0.25, m 3 ps ^Am^ is a good approximation. The neutrino mass spectrum 
can actually be determined to an acceptable degree of accuracy by using Eq. (13). 
For instance, we obtain m 3 m (3.8 — 6.1) x 10~ 2 eV, m 2 ~ (0.95 — 1.5) x 10~ 2 eV and 
mi ps (2.6 — 3.4) x 10 -3 eV, where x v ps 1/3 and y u ps 1/4 have typically been taken. 

3. Fig. 2 shows the outputs of sin 2 # 12 , sin 2 #23 an d sin 2 # 13 versus R v at the 3a level. 
One may observe that the maximal atmospheric neutrino mixing (i.e., sin 2 #23 ~ 0.5 
or sin 2 29 2 3 ~ 1) cannot be achieved from the isomeric lepton mass matrices under 
consideration. To be specific, sin 2 #23 < 0.40 (or sin 2 2# 2 3 < 0.96) holds in our ansatz. 
It is impossible to get a larger value of sin 2 #23 even if R u approaches its upper limit. 
In contrast, the output of sin 2 # 12 is favorable and has less dependence on R v . One 
may also see that only small values of sin 2 #13 (< 0.016) are favored. More precise 
experimental data on sin 2 623, sin 2 613 and R v will allow us to examine whether those 
parallel patterns of lepton mass matrices with six texture zeros can really survive the 
experimental test or not. 

4. Fig. 3 illustrates the results of two effective masses (m) e and (m) ee , three CP- violating 
phases (5, p, a), and the Jarlskog invariant J . It is obvious that (m) e ~ 10~ 2 eV for 
the tritium beta decay and (m) ee ~ 10~ 3 eV for the neutrinoless double beta decay. 
Both of them are too small to be experimentally accessible in the foreseeable future. 
We find that the maximal magnitude of J is close to 0.015 around 5 ~ 37r/4 (or 57r/4). 
As for the Majorana phases p and a, the relation (p — a) ps it/2 holds. This result 
is attributed to the fact that the matrix elements (a^b^c^) of U u are all imaginary 
and they give rise to an irremovable phase shift between V p \ and V P 2 (for p = e, p, r) 
elements through Eq. (23). Such a phase difference affects {m) ee , but it has nothing 
to do with (m) e and J . 

To relax the potential tension between the smallness of R v and the largeness of sin 2 d 2 3 , we 
shall incorporate a simple seesaw scenario in the six-zero textures of charged lepton and 
Dirac neutrino mass matrices in section IV. 
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B. Six Non-parallel Patterns (rank-3) 



The following six non-parallel patterns of M; and M l 
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(25) 



in which M v is of rank 3, are found to be compatible with current neutrino oscillation data 
at the 3a level. Given M\ and Mf being diagonalized respectively by the unitary matrices 
Ui and E 2 U U , where Ui jV = Pi, v 0^ v with Oi jV being simple functions of x l v and y l u as already 
shown in Eq. (23), the corresponding flavor mixing matrix reads 

V Pq = K)\V a + ( & p)* c ^ + , (26) 

where the subscripts p and q run respectively over (e, p, r) and (1, 2, 3), the phase parameters 
a and (5 are defined by a = (</?„ — ) — (20^ — 0;) and (5 = —{4> v + 4>i ), and an overall phase 
factor e*^ has been omitted. Taking account of the other five combinations of Mi and M v 
in Eq. (25), we notice that M/ n (for 1) and M[ n (for 2) can be diagonalized by E n Ui 
and (E n Ej)U u , respectively. Because the relation 

E T 2 {E X E T 2 ) = Ej(E 5 Ej) = Ej(E 6 Ej) = Ej(E 3 E?) = E^(E 4 E^) = E l (27) 

holds, Eq. (26) is universally valid for all six patterns. They are therefore isomeric. 

We do a numerical analysis of six non-parallel patterns of Mi and M v in Eq. (25). The 
parameter space of (x u ,y u ) or (a, (3) is found to be acceptable, when the 3a intervals of 
Afflj!, Airijj, sin 2 612, sin 2 ^23 and sin 2 #13 are used. Our explicit results are summarized in 
Figs. 4-6. Some brief discussions are in order. 

1. Fig. 4 shows the allowed parameter space of (x u , y v ) and (a, f3) at the 3a level. We see 
that (3 ~ (or j3 ~ 2ir) holds, while a is essentially unrestricted. Again, m 3 ps ^J~Krn^ 1 
is a good approximation. The neutrino mass spectrum can roughly be determined by 
using Eq. (13). Note that x v ~ 0.7 is marginally allowed - in this case, m\ and m 2 
are approximately of the same order. 

2. The outputs of sin 2 # 12 , sin 2 #23 and sin 2 6> 13 versus R v are illustrated in Fig. 5. We 
are unable to obtain the maximal atmospheric neutrino mixing (i.e., sin 2 #23 ~ 0.5 
or equivalently sin 2 2^23 ~ 1) from the non-parallel patterns of lepton mass matrices 
under consideration. Indeed, sin 2 #23 > 0.60 (or sin 2 2# 2 3 < 0.96) holds in our ansatz. 
It is impossible to get a larger value of sin 2 2^23 even if R v approaches its upper bound. 
In comparison, the output of sin 2 9 i2 is favorable and has less dependence on R v . Only 
small values of sin 2 13 (< 0.02) are allowed. 

3. The numerical results for (m) ee /m 3 versus (m) e /m 3 , J versus 5, and a versus p are 
shown in Fig. 6. Both (m) e ~ 10~ 2 eV and (m) ee ~ 10 -3 eV are too small to be 
observable. The maximal magnitude of J is close to 0.02 around 5 ~ ±vr/4, and the 
relation (a — p) ir/2 holds for two Majorana phases of CP violation. 

Comparing the parallel patterns of M^ u in Eq. (20) with those non-parallel patterns of M^ v 
in Eq. (25), we find that most of their phenomenological consequences are quite similar. 
Therefore, it is experimentally difficult to distinguish between them. 
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C. Twelve Non-parallel Patterns (mi = 0) 



Current neutrino oscillation data cannot exclude the possibility that the neutrino mass 
mi or m 3 vanishes. Hence M v is in principle allowed to take the rank-2 textures (M I? , M Ig 
and Mi g ) listed in Table II. After a careful analysis, we find that there exist four groups of 
non-parallel patterns of Mi and M v with m\ = 0, which are compatible with the present 
experimental data at the 3a level: 
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(28a) 

(28b) 
(28c) 
(28d) 



The possibility of m 3 = has been ruled out. With the help of Eqs. (18) and (19), it is 
easy to prove that three combinations of Mi and M v in each of the above four groups are 
isomeric. For the charged leptons, the expressions of (A h B h C{) and (a i; b i ,c i ) can be found 
in Eqs. (22) and (23). As for the neutrinos, we obtain 



A. 



m 3 (1 - Vv 



By = m 3 yfyy~ 



(29) 



where z v = Cy/m 3 . We see that it is impossible to fix C v (or B u ) in terms of m i: due to 
the fact that Det(M^) = holds. This freedom will be removed, however, once the flavor 
mixing parameters derived from Mi and My are confronted with the experimental data. To 
see this point more clearly, we write out the explicit results of nine elements of the lepton 
flavor mixing matrix V for every group of Mi and M u : 



V„ = K)Xe- + (b\)%e^ + (c l p yc u q 
with a = (fiy — (ip l — 4>i) and (3 = (p v — (p l corresponding to Eq. (28a); 

V m = (b l p )%e ia + {a\)%e^ + {c l p )% 
with a = (py — (pi and (3 = Cp v — ((pi — (pi) corresponding to Eq. (28b); 

V pq = {a l p Y~a^ + (c l p yb^ + {b l p )% 
with a = (by — ((fii — 20 ; ) and (3 = (p v + ; corresponding to Eq. (28c); and 



(30a) 



(30b) 



(30c) 



1 ~\*h"JP 



-q 



(30d) 
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with a = 4> v + 4>i and (3 = Cp^ — {<p l — 20 ; ) corresponding to Eq. (28d), where 



2fc - 2 2 



Vv ' 



i — ) u 2 ~ b I > u 3 ~ n — i ) 

3 = 0, c» = -i^ft=, % = -ft==- ( 31 ) 

V 1 + Vv V 1 + ^ 

In obtaining Eqs. (30c) and (30d), we have omitted an overall phase factor e~ % ^i . 

Note that the sum \a^\ 2 + \b^\ 2 (for % = 1,2,3) is independent of the free parameter z v . 
This result implies that V pq in Eq. (30b) can be arranged to amount to V pq in Eq. (30a). 

Indeed, the replacements z v •<=>- \jy v — z 2 and a (3 (or equivalently <fi u •<=>- (p v ) allow 
us to transform (V pl: V p2: V p3 ) of Eq. (30a) into (— V pl: V p2 , V p3 ) of Eq. (30b). The extra 
minus sign of V p \ appearing in such a transformation does not make any physical sense, 
because it can be removed by redefining the phases of three charged lepton fields. Thus we 
expect that Eqs. (30a) and (30b) lead to identical results for lepton flavor mixing and CP 
violation. One may show that Eqs. (30c) and (30d) result in the same lepton flavor mixing 
and CP violation in a similar way. For this reason, it is only needed to numerically analyze 
the non-parallel patterns of Mi and M v in Eqs. (28a) and (28c). 

A numerical analysis indicates that the parameter space of (y u , z„) or (a, (3) can be found, 
if the 3<t intervals of Am^, Aro^, sin 2 9u, sin 2 6*23 and sin 2 #13 are taken into account. Our 
results are summarized in Figs. 7-12. Some comments are in order. 

1. The parameter space and predictions of Mi and M v listed in Eq. (28a) are shown in 
Figs. 7-9. We see that (3 ~ re is favored but a ~ tc is disfavored. The neutrino mass 
spectrum has a clear hierarchy: x v = and y u ~ 0.25. The outputs of sin 2 6> 12 and 
sin 2 023 are well constrained, and they seem to favor the corresponding experimental 
lower bounds. Again, it is impossible to obtain the maximal atmospheric neutrino 
mixing. We observe that large values of sin 2 #13, more or less close to its experimental 
upper limit, are strongly favored. This interesting feature makes the present ansatz 
experimentally distinguishable from those given in Eqs. (20) and (25). As a straight- 
forward consequence of the normal neutrino mass hierarchy, the results of (m) e and 
(m) ee are both too small to be observable in the near future. The maximal magnitude 
of J is close to 0.02 around \b\ ~ ±7r/7. As for the Majorana phases, we get the 
relation (a — p) it/2 (or — 3ir/2). 

2. The parameter space of Mi and M v in Eq. (28b) can be obtained from Fig. 7 with 



the replacements z v •<=>- yy u ~ z l and a •<=>- f3. Such replacements are actually 

equivalent to B v -<=>- C v and <p v •<=>- <p v between M v in Eq. (28a) and its counterpart 
in Eq. (28b). The phenomenological consequences of Mi and M u in both cases are 
identical, as already shown above. 
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3. Figs. 10-12 show the allowed parameter space and predictions of Mi and M v listed in 
Eq. (28c). We see that a ~ n and (3 ~ (or 2n) are essentially favored. The neutrino 
mass hierarchy is quite similar to that illustrated in Fig. 7. The output of sin 2 #23 seems 
to favor the corresponding experimental upper bound, and the maximal atmospheric 
neutrino mixing cannot be achieved. In comparison, the outputs of sin 2 6*12 and sin 2 9 13 
are favorable and have less dependence on R v . Note that the predictions of this ansatz 
for (m) ee and J may reach OAm^ (at (m) e ~ 0.15m3) and 0.03 (at 5 ~ ±37r/4), 
respectively. Both results are apparently larger than those obtained above. Again, the 
relation (a — p) ps it/2 (or — 37r/2) holds for two Majorana phases. 

4. The parameter space of Mi and M v in Eq. (28d) can be obtained from Fig. 10 with the 
replacements z v •<=>- \Jy u — z"l and a •<=>- (3. Their phenomenological consequences 
are identical to those derived from Mi and M v in Eq. (28c). 

The main unsatisfactory output of twelve non-parallel patterns of Mi and M u , just like the 
one of six parallel patterns of M; and M v in Eq. (20), is that sin 2 2# 2 3 cannot reach the 
experimentally-favored maximal value. Whether this is really a problem remains to be seen, 
especially after more accurate neutrino oscillation data are accumulated in the near future. 



IV. A SEESAW ANSATZ OF LEPTON MASS MATRICES 

To illustrate, let us discuss a simple way to avoid the potential tension between the 
smallness of R u and the largeness of sin 2 #23 arising from those parallel patterns of Mi and 
M v in Eq. (20). In this connection, we take account of the Fukugita-Tanimoto-Yanagida 
hypothesis [23] together with the seesaw mechanism [24] - namely, the charged lepton mass 
matrix M t and the Dirac neutrino mass matrix M D may take one of the six parallel patterns, 
while the right-handed Majorana neutrino mass matrix Mr takes the form Mr = M Ei with 
Mq denoting a very large mass scale and E\ being the unity matrix given in Eq. (18). Then 
the effective (left-handed) neutrino mass matrix M v reads as 

M 2 

M v = MvM^Ml = -f- . (32) 

Mq 

For simplicity, we further assume M D to be real (i.e., <p D = ip D = 0) . It turns out that the 
real orthogonal transformation Ud, which is defined to diagonalize M D , can simultaneously 
diagonalize M v : 

UEM„U D = {UlM *y 2 = I m 2 , (33) 

where mj = df/M with di standing for the eigenvalues of M D . In terms of the neutrino 
mass ratios x v = mi/m 2 = (rfi/rf 2 ) 2 and y u = m 2 /m3 = (c^/cfa) 2 , we obtain the explicit 
expressions of nine matrix elements of U u = Un- 
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1/2 



On = - 



+ 



&? = + 



'i + v/^)(i - v /x~^;)(i - ^ + v /s^;) 

v /^J(l + yjx v y v ) 



(1 + + v^)(l - v/^J + y/x v y v ) 

y vy /x v y v {l - y/x^) 

.(1 - y/x v y v ){l + y/y^)(l - yfWv + y/x v y v ) 

y/xz{i-y/y7) 



1/2 



(1 + y/x^)(l - ^x u y u ) 



= + 



= + 



1 + yjx v y u 



[i + y/xj(i + v^;) 



1/2 



1/2 



Ci = - 



C U 2 



(1 - v/x^)(l + v /y^)_ 

y/x v y v {l - y/xZ)(l + sjx v y v ) 



1 + 



v /^;(i - v /x^)(i - 



1/2 



<y. 



y v + yjx v y v 

nl/2 



+ 



(1 + v /x^)(l + v/^J)(l - v^J+ ^/x u y u )_ 

(1 - y/y^Kl + y/gj/jJ) 

'1 - y/x v y v ){\ + v/^;)(l - + y/x v y v ) 



1/2 



(34) 



The lepton flavor mixing matrix V = U}U U remains to take the same form as Eq. (24), but 
the relevant phase parameters are now defined as a = — tp t — (3 and f3 = — <p l . Comparing 
between Eqs. (23) and (34), one can immediately find that the magnitudes of (9 12, #23, #13) 
in the non-seesaw case can be reproduced in the seesaw case with much smaller values of 
x v and y v . The latter will allow R v to be more strongly suppressed. It is therefore possible 
to relax the tension between the smallness of R v and the largeness of sin $23 appearing in 
the non- seesaw case. A careful numerical analysis of six seesaw-modified patterns of lepton 
mass matrices does support this observation. The results of our calculations are summarized 
as follows. 

1. We find that the new ansatz are compatible very well with current neutrino oscillation 
data, even if the 2a intervals of Am^, Ara^, sin 2 12 , sin 2 #23 and sin 2 13 are taken into 
account. Hence it is unnecessary to do a similar analysis at the 3a level. The parameter 
space of (x u ,y u ) and (cx,f3) is illustrated in Fig. 13, where x v ~ y v ~ 0.2 and (3 ~ n 
hold approximately. Ag approximation. The values of three 
neutrino masses read explicitly as m 3 pa (4.2 — 5.8) x 10~ 2 eV, rri2 ~ (0.84— 1.2) x 10~ 2 
eV and m 1 ~ (1.6 — 1.9) x 10 -3 eV, which are obtained by taking x v ~ y v ~ 0.2. 

2. The outputs of sin 2 # 12 , sin 2 #23 and sin 2 # 13 versus R v are shown in Fig. 14 at the 2a 
level. One can see that the magnitude of sin 2 12 is essentially unconstrained. Now 
the maximal atmospheric neutrino mixing (i.e., sin 2 #23 ~ 0.5 or sin 2 29 23 ~ 1) is 
achievable in the region of R u ~ 0.036 — 0.047. It is also possible to obtain sin 2 9 13 < 
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0.035, just below the experimental upper bound [4]. If sin 2 26 13 > 0.02 really holds, 
the measurement of #13 should be realizable in a future reactor neutrino oscillation 
experiment [25]. 

3. Fig. 15 illustrates the numerical results of (m) e , (m) ee , 5, p, a and J . We obtain 
(m) e ~ 10~ 2 eV for the tritium beta decay and (m) ee ~ 10~ 3 eV for the neutrinoless 
double beta decay - both of them are too small to be experimentally accessible in the 
near future. We see that \ J\ ~ 0.025 can be obtained. Such a size of CP violation is 
expected to be measured in the future long-baseline neutrino oscillation experiments. 
As for the Majorana phases p and a, the relation a ~ p holds. This result is easily 
understandable, because U v is real in the seesaw case. It is worth mentioning that the 
effective neutrino mass matrix M v does not persist in the simple texture as M\ has, 
thus the allowed ranges of 8, p and a become smaller in the seesaw case than in the 
non-seesaw case. 

It should be noted that the eigenvalues of Md and the heavy Majorana mass scale M 
are not specified in the above analysis. But one may obtain \di/d 2 \ = ^/x^ ~ 0.4 and 
l^/cfal = \[Vv ~ 0-4- Such a weak hierarchy of \d 2 \, |d 3 |) means that M D cannot 
directly be connected to the charged lepton mass matrix M ; , nor can it be related to the 
up-type quark mass matrix (M u ) or its down-type counterpart (M<j) in a simple way. If the 
hypothesis M R = MqEi is rejected but the result UjM u U u = Diag{m l5 m 2 , %} with U u 
given by Eq. (34) is maintained, it will be possible to determine the pattern of M R by means 
of the inverted seesaw formula M R = M-^M^Md [26] and by assuming a specific relation 
between Md and M u . For example, one may simply assume Md = M u with M u taking the 
approximate Fritzsch form, 



M u 



/ yjm u m c \ 

^m u m c ^m c m t 
\ ^m c m t m t J 



(35) 



Just for the purpose of illustration, we typically input x u ~ y u ~ 0.18 as well as m u /m c 
m c /m t ~ 0.0031 and m t « 175 GeV at the electroweak scale [6]. Then we arrive at 



M R ~ 3.0 x 10 15 x 



/6.1xl0~ 8 1.2 x 10~ 5 2.0xl0~ 4 \ 

1.2 x 10~ 5 3.5 x 10~ 3 5.9 x 10^ 2 
\2.0xl0- 4 5.9 x 10~ 2 1 ) 



(36) 



in unit of GeV. This order-of-magnitude estimate shows that the scale of M R is close to that 
of grand unified theories A GUT ~ 10 16 GeV, but the texture of M R and that of M D (or M t ) 
have little similarity. It is certainly a very nontrivial task to combine the seesaw mechanism 
and those phenomenologically-favored patterns of lepton mass matrices. In this sense, the 
simple scenarios discussed in Refs. [22,23] and in the present paper may serve as a helpful 
example to give readers a ball-park feeling of the problem itself and possible solutions to it. 

Of course, a similar application of the seesaw mechanism to the non-parallel patterns of 
lepton mass matrices is straightforward. In this case, an enhancement of sin 2 26* 2 3 up to its 
maximal value can also be achieved. 
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V. SUMMARY 



To summarize, we have analyzed 400 combinations of the charged lepton and neutrino 
mass matrices with six texture zeros in a systematic way. Only 24 of them, including 6 
parallel patterns and 18 non-parallel patterns, are found to be compatible with current neu- 
trino oscillation data at the 3a level. Those viable patterns of lepton mass matrices can be 
classified into a few distinct categories. The textures in each category are demonstrated to 
have the same phenomenological consequences, such as the normal neutrino mass hierarchy 
and the bi-large flavor mixing pattern. We have also discussed a very simple way to incor- 
porate the seesaw mechanism in the charged lepton and Dirac neutrino mass matrices with 
six texture zeros. We illustrate that there is no problem to fit current experimental data 
even at the 2a level in the seesaw case. In particular, the maximal atmospheric neutrino 
mixing can naturally be reconciled with a relatively strong neutrino mass hierarchy. Our 
results for effective masses of the tritium beta decay and the neutrinoless double beta decay 
are too small to be experimentally accessible in both the seesaw and non-seesaw cases, but 
the strength of CP violation can reach the percent level and might be detectable in the 
upcoming long-baseline neutrino oscillation experiments. 

We conclude that the peculiar feature of isomeric lepton mass matrices with six texture 
zeros is very suggestive for model building. We therefore look forward to seeing whether 
such simple phenomenological ansatze can survive the more stringent experimental test or 
not in the near future. 

This work was supported in part by the National Natural Science Foundation of China. 
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TABLES 



TABLE I. The best-fit values, 2a and 3a intervals of Am^, lAmJ-J, sin 2 #i 2 , sin 2 #23 and 
sin 2 obtained from a global analysis of the latest solar, atmospheric, reactor and accelerator 
neutrino oscillation data [14]. 





Am^ (10- 5 eV 2 ) 


\Am 2 3l \ (10- 3 eV 2 ) 


sin 2 612 


sin 2 #23 


sin 2 9 n 


Best fit 


6.9 


2.6 


0.30 


0.52 


0.006 


2a 


6.0-8.4 


1.8-3.3 


0.25-0.36 


0.36-0.67 


< 0.035 


3a 


5.4-9.5 


1.4-3.7 


0.23-0.39 


0.31-0.72 


< 0.054 



TABLE II. The type-I textures of a symmetric lepton mass matrix M (i.e., Mi or M u ) and the 
corresponding forms of the phase matrix P (i.e., Pi or P v ) and the unitary matrix O (i.e., Oi or 
O v ) used to diagonalize M, in which (^4, B, C) or (A, B, C) are defined to be real and positive. 
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FIGURES 



Q.2 D.3 Q.4 OS 



FIG. 1. Parallel patterns of M/ and M v in Eq. (20): the parameter space of (x u ,y u ) and (a, 
at the 3cr level. 
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FIG. 2. Parallel patterns of M/ and M v in Eq. (20): the outputs of sin 2 12 , sin 2 6> 2 3 and sin 2 13 
versus R v at the 3a level. 



23 




FIG. 3. Parallel patterns of Mi and M v in Eq. (20): the outputs of ((m) e , (m) ee ), (5, J) and 
(p, a) at the 3a level. 
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FIG. 4. Non-parallel patterns of M/ and M v in Eq. (25): the parameter space of (x u ,y u ) and 
(a, /3) at the 3<r level. 
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FIG. 5. Non-parallel patterns of Mi and M;, in Eq. (25): the outputs of sin 2 8u, sin 2 823 and 
sin 2 versus R v at the 3d level. 
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FIG. 6. Non-parallel patterns of M; and M v in Eq. (25): the outputs of ((m) e , (m) ee ), (<5, J") 
and (p, <t) at the 3<r level. 
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FIG. 7. Non-parallel patterns of Mi and M v in Eq. (28a): the parameter space of (y v ,z v ) and 
(a, /3) at the 3a level. 
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FIG. 8. Non-parallel patterns of M; and M v in Eq. (28a): the outputs of sin 2 6\2, sin 2 #23 and 
sin 2 versus R v at the 3a level. 
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FIG. 9. Non-parallel patterns of Mi and in Eq. (28a): the outputs of ((m) e , (m) ee ), (5, J") 
and (p, a) at the 3<r level. 
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FIG. 10. Non-parallel patterns of Mi and M v in Eq. (28c): the parameter space of (y u , z v ) and 
(a, (3) at the 3a level. 
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FIG. 11. Non-parallel patterns of Mi and M v in Eq. (28c): the outputs of sin 2 612, sin 2 #23 and 
sin 2 versus R v at the 3<r level. 
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FIG. 12. Non-parallel patterns of M; and M v in Eq. (28c): the outputs of ((m) e , (m) ee ), (5, J") 
and (p, a) at the 3<r level. 
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FIG. 13. A simple seesaw example: the parameter space of {x v ,y v ) and (a, (3) at the 2a level. 
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FIG. 14. A simple seesaw example: the outputs of sin 2 #12, sin 2 023 and sin 2 #13 versus R v at 
the 2(7 level. 
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FIG. 15. A simple seesaw example: the outputs of ((m) e , (m) ee ), (5, J) and (p,a) at the 2a level. 
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